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I. INTRODUCTION 

Electromagnetic form factors of hadrons contain important information about the intrin- 
sic structure of these particles and about their coupling to the external electromagnetic field. 
This electromagnetic structure has been explored mainly by electron scattering experiments, 
giving information on the electromagnetic form factor in the space like momentum region 
and thus on the spatial distribution of charges and magnetic moments inside the hadrons 
in general, and the nucleon in particular. The available data cover the momentum range 
from q 2 = up to about q 2 = 35 GeV 2 and can be well parametrized by a dipole fit [|TJ. 
There exist models who explain the data by assuming a certain number of poles which all, 
unfortunately, are in the so-called "unphysical" region (not accessible by experiments on 
on-shell nucleons) 0||. 

The data are much more sparse for time like momentum transfers, where the excitation 
of the nucleon and its decay is studied. Here the only available data come from experiments 
at LEAR, exploiting pp annihilation |4]]. These data thus naturally start at momentum 
transfers larger than two times the mass of the proton. They show that conventional pole 
fits e. g. cannot be applied in this region any more. It is thus of special interest, to 
explore the "unphysical region" and see whether there exists a rich structure of poles and 
thus whether the so-called Vector Meson Dominance (VMD ||), which is underlying these 
pole fits, is a universal property of all hadrons or if it holds only for pions. 

The region around momentum transfers corresponding to the vector meson mass (around 
750 - 800 MeV for p and u mesons) is so interesting because these pole fits predict here a 
very pronounced resonance structure in the form factors ||. Since access to this region 
for on-shell protons is impossible, the only alternative is to look for half off-shell processes. 
Indeed, dilepton production in hadronic collisions (bremsstrahlungs-dileptons) offers access 
to the half off-shell electromagnetic form factor in the time like region || . 

Experiments of this kind have been performed by the DLS collaboration for p + p, p + A 
and A + A collisions M. Simulations of these processes show for p + p and p + A a clear 



window for these bremsstrahlungs contributions where these are not overshadowed by other 
processes ||; for A + A a strong u peak is predicted 0. 

Under the assumption of VMD the elementary dilepton production processes NN — > 
NNe + e~, tcN — > Ne + e~ and jN — ► Ne + e~ were studied by the Giessen group [p],[10|,[TT 



It is, however, unclear if the time like electromagnetic form factor is influenced by the off- 
shellness of the intermediate nucleon. The purpose of this paper is to model the relevant 
vertex, to study its off-shell dependence and to investigate if VMD is still visible if this 
off-shell dependence is properly taken into account. The calculations are thus meant to 
stimulate and to provide some guidance for experimental investigations of this important 
hadronic property. 

To obtain the general vertex, one needs a dynamical model that describes the electro- 
magnetic structure of the nucleon. It would be most desirable to use the quark degrees of 
freedom. Unfortunately, state of the art quark models of the nucleon do not allow to study 
off-shell effects and excitations quantitatively, but only in a qualitative way. Their success 
is mainly confined to space like properties of real nucleons at the present time |H| . In this 



paper, therefore, all calculations are performed using hadronic degrees of freedom, follow- 
ing the concept of Naus and Koch |T3| or Tiemeijer and Tjon [TJJ], who performed similar 
calculations for the space like regime. 



From a spectral analysis of form factors [15] one has learned about two important in- 



gredients: 7riV-scattering and 7T7r-scattering. While 7T7r-scattering is resonant, and thus in 
this channel best described by taking into account a coupling to the p-meson, 7riV-scattering 
is more difficult to implement, and will be handled only in a very schematic way in the 



present model. On the other hand, semiphenomenological models [|T6| , |T7|1 show the success 
of a description of the form factors in a cloud/core picture. This will be discussed later in 
this paper. 

Constrained by gauge invariance of the electromagnetic interaction, the vertex is con- 
structed from a pion loop expansion of the nucleon propagator by coupling external photons 
to each charged particle in the loop. In this loop expansion the A-resonance is taken into 



account because, first of all, it is necessary for the reproduction of 7riV-scattering, and sec- 
ondly, it is the first important resonance to contribute to the off-shell effects. In addition 
to that, from an analysis of the Gerasimov-Drell-Hearn sum rule |[L8| , the A is found to be 
most important to understand the anomalous magnetic moments of the nucleons. 

This paper is organized as follows. In section [TI], the general structure of form factors 
for nucleons and pions is discussed, gauge invariance and its consequences are introduced. 



In section JTJ, the model is presented, and gauge invariance is proven for a simple case. A 



discussion of VMD in this context is given in section |V|. Results and comparison to data 



are shown in section [V| Section [Vl| concludes the paper. 



II. GENERAL STRUCTURE OF FORM FACTORS 

A. Form Factor of Nucleons 

The most general form of the electromagnetic interaction vertex for nucleons can be split 
into an isoscalar part and an isovector part 

Y^p\p)=Y^p',p) + r^ v {p',p) . (1) 

The vertices of proton or neutron are then linear combinations T^(p',p) = Tg(p',p) + 

iW,p), W, P ) = r$(jpf, P ) - iw,p). 

The isoscalar or isovector parts can be split up further. Let p denote the 4-momentum 
of the incoming nucleon, p' the 4-momentum of the outgoing nucleon, and q = p — p' the 
4-momentum of the outgoing photon. This choice is more convenient for the case of time like 
momentum transfer and for dilepton production, but is different to what is conventionally 
used in the literature, where q is the 4-momentum of the incoming photon; this amounts to 



some sign changes in the decomposition of the vertex compared to other papers ||13| , [14| . [19 



According to |T9| the isoscalar/isovector vertex can be expressed as 



n, v (p',p) = e E {7 • v'Y (Ki v r - Klv^- - *£&<?) b • vY ■ (2) 

r',r=0 



The 24 functions A-g V (i = 1..3; r, r' = 0, 1) are scalar functions of the three variables p 2 , p' 
and q 2 . Introducing a shorthand notation 

0? = j>*, 0% = — —, O^ = -q", (m = nucleon mass) 

and the projection operator 

a s (p) = S1 'l^ - withiy = v /^, s = ±l , (3) 

(0) can be written as (drop S and V for convenience) 

PV,P)=e£ £ A s V)^(^',^;g 2 )0fA*(p) . (4) 

i=l s,s'=±l 

The functions Ff s are linear combinations of the A\ r . 

Since for p 2 = m 2 the operator A s (p) is just the projection operator to positive (s = +1) 
or negative (s = —1) energies, r M (p',p) taken between on-shell spinors reduces to 
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u(p')T >1 (p',p)u(p) = ey^ j u{p')F^ r+ [rn^rn] q 2 )Q^u{p) . 

i=l 

It turns out that F 3 ++ (m, m; q 2 ) = (time reversal invariance), so the on-shell vertex takes 
on the well known form |2"0|| 

r ;„uv„ -\ 

u(p) . 



u(p')T»(p',p)u(p) = eu(p') F+ + {q 2 )Y - F 2 ++ (q 2 ) ' 



2m 

At q 2 = the form factor -F\ ++ measures the electric charge and -F 2 ++ measures the anomalous 
magnetic moment. 

The vertex T fl (p',p) contains self energy corrections in the external legs. For this reason 
it is often called reducible vertex. In this paper it is also called full vertex in order to indicate 
that in experimental measurements these self energy corrections are always included. To 
eliminate these corrections, the irreducible vertex is defined by 

S(p'M. r (p',p)S(p) = S (p')T»(p' 1 p)S (p) . (5) 

For the irreducible vertex a decomposition exists similar to the one for the full vertex, as 
will be proven now. In @ S (p) is the free propagator, S(p) is the full propagator, including 



all self energy corrections E(p) = Sq 1 (p) — S^ 1 (p). Most generally E(p) is decomposed into 
a vector part Ey(W / ) and a scalar part Eg(W) by E(p) = Ey(VT)7 -p — E 5 (W / )m, so the 
full propagator takes on the form 

S-\ P ) = 7 ■ p - m - E(p) = 7 • p(l - E V {W)) - m(l - E S {W)) = 7 • p* - m* . 

It proves helpful to introduce the positive and negative energy projections of the self energy 
and the propagator E ± (W / ) and 5 ,± (W / ) by 

E(p) = £ A s (p)Z s (W), S(p) = £ A*(p)S s (W) . 

s=±l s=±l 

Using the properties of the projection operators 

A + (p) + A-(p) = l, A+(p)-A-(p) = l^ , (6) 

one finds for the self energy 

Y, S (W) = sY, v (W)W -^ s (W)m (7) 

and for the full propagator 

1 1 

( > ~ sW * _ m * ~ sW _ ~ _ y, s (W) ' ^ ' 

where W* = W(l — Ey(W)) and m* = m(l — Hs(W)). With this notation one finds further 

s Q (p) = Y,A s (p)s s (wy, s°(w) = — i — (9) 

g sW — m 

S" 1 ^) = £A s (p)5- ls (lU); S" 1 ^^) = sW* - m* (10) 

s 

S 1 (p)=J2A s (p)S ls (W) ] S ls (W)=sW-m . (11) 

s 

Now one can easily derive the relation between the reducible vertex (^) and the irreducible 
vertex (H): 

rf rr ( P » = s- 1 (p')s (p')r'(p',p)s (p)s- 1 (p) 

*=1 s's 

= eEE AS '(p')//' s (^',W;g 2 )OrA s (p) . (12) 

8=1 s's 



Eq. fljjD represents the desired decomposition of the irreducible vertex in terms of the irre- 
ducible form factors ff s (W, W; q 2 ). 

The Ward-Takahashi-identity (WTI) relates the irreducible vertex to the full propagator 



m- 



q,K r (p',P) = eQ(S- 1 (p)-S- 1 (p')) . (13) 

Q is the appropriate charge operator. 

Two useful identities will be introduced here: 

A S V)7 • qA s {p) = A s '{p'){sW - s'W')A s {p) (14) 

-A s '(p')ia^A s (p) = A s '(p')[(sW + s'W') lti - (p + p') ,]A S (p) . (15) 



(|I~5|) is a general form of the Gordon-identity J20f . With the use of ( |T4| ) and q^O^ = the 
projection of (|T3| ) onto positive and negative states yields 

f('{sW - s'W) - q 2 f°' s = Q (S- 1S (W) - S- 1S '(W')] 

= Q ((sW* - s'W'*) - (m* - m'*)) . (16) 



For the case of an outgoing on-shell particle (W = m) of positive energy (s' = +1) this 
reduces to 

sW — m sW — m 



or for the full form factor y (|17|) • f^_ 



a 2 

F+ s F+ s = Q . (18) 

1 sW* -m* 3 K J 

To get flT8| ) one needs further ^*Zm>* , = 1, i. e. the full propagator has a pole of unit 
residue at the physical mass m of the nucleon. 

For real photons (q 2 = 0) one recovers that F^ s = 1 for protons and F^ s = for 
neutrons. This holds for the full form factor even if the incoming nucleon is off-shell. Note 
that fllTD implies that fi~ s (q 2 = 0) depends on the nucleon's off-shellness and only reduces to 



the real charge for on-shell incoming particles. The latter must be true because for on-shell 
particles there is no difference between full and irreducible vertex. 

Note that the WTI does not pose any constraint on the magnetic form factors F^~ s . Note 
also, that F^~ s (m,m]q 2 ) = f^ s (m,m]q 2 ) have to vanish for all q 2 in order to obtain finite 
contributions on the lhs of ( |T7D and flT8|). 

To obtain the full vertex, the knowledge of the full propagator is needed, which can be 
obtained from /{*" s using relation ( |17|) taken at q 2 = 0: 

S- 1S (W) = sW* - m* = /ij(m, W; q 2 = 0)(sW - m) . 

The index p stands for proton. This allows to write the half off-shell full form factor as 

It is thus sufficient to calculate the irreducible vertex only. 



From (16) more relations can be derived especially for q — 0: 

f+ + (W,m ] 0) = f+ + (m,W ] 0) 
f+-(m,W-0) = f 1 -+(W,m;0) 
fx~{m, m; 0) = 1, f{~(m, m; 0) = 

So the off-shell full form factor is 

Ffw wq 2 ) = , agjg 



(20) 



Experiments always measure the full form factor. F 3 S s is never accessible by experiments 
since O^j^ = for any conserved current j^. 



B. Form Factor of Pions 



The most general form of the pion photon vertex is [22 



W,p) = eQMW ,P 2 ; Q 2 )P >fl + W\V 2 \ q 2 )p 



A more convenient and also better known notation is 

with P L =p' + p, P T = P L - q{P L ■ q/q 2 ), q = p - p' . 
If jn is a conserved current, one measures 

j,K(p',p) = eQ n [A(p' 2 ,p 2 ; q 2 ) + B(p' 2 ,p 2 ; q 2 )} 3 ,P£ , 

i. e. the sum of A and B. The WTI requires 

q,K(p',p) = eQ^P£A{p>\p 2 - q 2 ) = eQ.^V) - D^(p' 2 )} . 

Thus only A is constrained by the WTI. 

The on-shell form factor of pions is a measurable quantity. It is given to good precision 
by the vector meson dominance (VMD) hypothesis [p3fl . In fact the VMD works so well 



that one is led to assume that the bare pion is essentially a structureless particle p4 |. It 



is therefore safe to neglect effects of the pion's off-shellness and to assume that A and B 
depend on q 2 only: 

J^(p',p) = eQ 7r F 7T (q 2 ) 3lM P^ , 

where F n {q 2 ) is the measured form factor. This fixes B(q 2 ): 

B{q 2 ) = F^{q 2 ) - A{q 2 ) . 

Given any form of the pion self energy, this allows to maintain gauge invariance as well as 



the measured form factor. Section [TV] will explain the VMD hypothesis and will show how 
to carry it over to the nucleon. 

III. MODEL FOR THE FORM FACTOR 

As mentioned in the introduction, there exists a subtle interplay of tcN- and tctt scat- 
tering in describing the electromagnetic properties on nucleons. This section is devoted to 
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introduce the part of the model that is suited to describe the irN interaction. Subsection 
1 11 A| introduces the Lagrangians for nucleons, pions and As. In a naive picture the coupling 



of photons is introduced. Subsection |111 B| describes the coupling of the photon to these 



fields in the correct way and shows that this is equivalent to the picture developed in sub- 
section |111 A| , which is thus respecting the local U(l)-symmetry of QED. In section [IV] the 



7T7T interaction is modeled in terms of vector meson dominance. 

A. Interaction of Mesons, Baryons and Photons 

Part of the structure of the nucleon is due to the meson cloud that dresses the bare 
nucleon. This idea is well established and several phenomenological models exist that take 
it into account [16,17||. For electromagnetic form factors Naus et al. |13[], Tiemeijer et al. 



T4J or Bos et al. |25[ performed detailed calculations based on a meson baryon interaction 
picture, all giving essentially the same results for the shape and the off-shell dependence of 
the form factors, but differing somewhat in their predictions for the anomalous magnetic 
moments. However these authors restricted themselves to space like momentum transfer and 
to off-shell nucleons with W < m + m n in order to avoid poles due to decay into inelastic 
channels. Since the purpose of this paper is to compute the form factor in the time like 
region, the decay modes have to be included; also the A-resonance will be considered in 
the calculation because for off-shell nucleons with W > m the A is no longer kinematically 
supressed. 

The model used here is based on the Lagrangian densitiy for pions and nucleons with 
pseudovector coupling. 

C N = f (7 ■ p - m)tf + \[{d^)\d^) - m^vf] + ^g^ 7 57 M r^ . (21) 



Including the A leads to an additional term in fl2T|): 



C A = ^ A A^l + ^ («r£f¥0 M 7? + h.c.) (22) 



2m 
with 
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A^ = (l-p-m A )g^ v - {l^Pu+P^lu) +1^1 -piu + mAlvlu , (23) 



as derived in the Rarita-Schwinger formalism fl26|j2Tf , and T being the matrix that couples 



isospin 3/2 to isospin 1/2 © 1. The A is here treated as a stable particle; its finite decay 
width must be neglected at the order of diagrams discussed here. Using a momentum 
dependent width would amount to including self energy corrections to the A propagator; 
this corresponds to diagrams of higher order in pion lines. Also a whole set of new diagrams 
would be necessary to maintain gauge invariance. On the other hand a constant width in 
the A propagator always yields complex form factors, even for on-shell nucleons, because 
the relevant thresholds are not taken into account. 

This subsection will describe how the corrections to the electromagnetic interaction ver- 
tex can be constructed for diagrams including one pion loop. Starting point will be the 
nucleon propagator, which is up to the one pion loop level given by fig. |l[ In the naive 
picture the photon couples to all charged particles individually. That this indeed fulfills the 
WTI will be proven in the following subsection. The irreducible vertex is given by fig. |2|, 
where the coupling of the photon to the hadrons is according to the usual Feynman rules. 
Since the pion coupling is chosen to be pseudovector, additional contact terms arise (fig. 
|d,e). 

Since all the loop diagrams diverge, a regular ization by a covariant cutoff of the form 

fM 2 ) = "f^ = (ml - tf)D c {q) (24) 

is used. This form of the cutoff is chosen because it can be visualized as the propagator 
D c (q) of a particle of mass A with the same quantum numbers as the pion; this ficticiuos 
cutoff-particle will be called " cop" . The iViW- vertex is regulated with the monopole form, 
whereas the iVA7r- vertex needs a dipole cutoff to yield convergence because the A-propagator 
is 0(p). In the spirit of the previous paragraph the photon will then also couple to the cop if 
the intermediate pion carries charge. The photon-pion coupling must therefore be replaced 
by a sum of vertices given in fig. Q the double dashed line denotes the propagator of the 
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cop. The diagrams with an internal A-propagator have two cop propagators because of the 
dipole cutoff. The pion-photon-vertex in this case must thus be replaced by five diagrams. 
The mass of the cop must be chosen large enough so that unphysical decay modes are 
avoided. The possible decay modes and thus contributions to the imaginary part of the 
form factor can be found by applying the usual cutting rules. If the incoming nucleon is 
far enough off-shell (W > m + m n ) it can decay into an on-shell nucleon and a real pion 
as depicted in fig. f|a. In principle the off-shell nucleon can also decay into a nucleon and a 
cop if W > m + A. This inequality sets the limit of the model. For the case of this paper 
W max will be restricted to 2 GeV, so A should be somewhat larger than 1 GeV. This is in 
good agreement with cutoff values used in meson exchange potentials. If the cut is taken as 
indicated in fig. |]b it becomes clear that there is also a restriction on the invariant mass M of 
the photon. If M = Jqi > 2m !n 7r + 7r~-production becomes possible. Pion-cop-production 
is not possible if M is restricted to M < m n + A ps m^ + W max -mwl GeV. 

B. Ward-Takahashi-Identity 

This subsection is devoted to the proof of correctness of the picture developed above. 
First, the appropriate electromagnetic interaction vertex will be derived from the La- 



grangians (pl , p2|) by minimal coupling. Then the "reduced" formalism will be introduced 
to keep track of the cop. Finally it will be proven that the WTI is fulfilled. 

Since the photon is the gauge field of local U(l) symmetry, it is introduced into the 
Lagrangian by substituting p^ — ► p^ — eQA^, where Q is the charge operator. This leads to 
the usual interaction terms for nucleons and pions: 

to a contact term because of the derivative in the pseudovector coupling of pions and nucle- 
ons: 



9NN-K . 

2m 



C = -e^VTeffVQ^Ap , 
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and to terms from £&: 



gNAn 



—e- 



2m 



(# A f ^4^ + h.c. 



The operators Qn,N,A return the appropriate charge. Since the A can decay into a 7 and a 
nucleon, an additional term must be included in the Lagrangian: 

£ A ->n 7 = i^(* A7 57^ 3 V|/i^ + h.C.) . 



The corresponding vertex is 75(7'^ — 7 • qg fJ,u )Gi where the factor a/2/3 from T 3 and 



the coupling constant have been absorbed in G\. According to |28| this vertex is mainly 
responsible for the Ml multipole which dominates the decay, q is the momentum of the 
outgoing photon. Note that the index /1 of the above vertex contracts with the A, while the 
index v contracts with the photon field. 

To prove gauge invariance of the model developed in this paper, three steps need to be 
done. To obtain a gauge invariant coupling to the pion, the coupling to the cop needs to be 
investigated, as well. So, first, the cop-photon interaction vertex is constructed, then, as a 
second step, the effective pion-photon vertex is constructed to fulfill the WTI locally using 



the "reduced" formalism of Gross and Riska |P2 |. The third step is to actually prove the 
WTI for the photon-nucleon interaction. 

Since the cop is not a fundamental particle, it does not appear in the Lagrangian. There- 
fore, so far it does not interact with the electromagnetic field. To obtain gauge invariance 
it needs the same interaction vertex as the pion: 

T^(p\p) = T^(p',p) = eQ^p + p'Y , 

where p and p' are the momenta of the incoming and outgoing particles respectively. Note 
that the so constructed vertex obviously fulfills the WTI of a free particle: 

(P - p'),T^(p',p) = eQ^P 2 ~ P' 2 ) = eQ* {Dc\p) ~ £> C ~V)) • (25) 
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Next the "reduced" formalism is introduced. This part of the paper follows closely fl22 



who treat the cutoff function f c as a contribution to the pion polarization function, which is 
possible because f c is a function of the pion momentum only. The "reduced" pion propagator 
D R (p) is defined as 

d r (p) = fc(p 2 )DMf c ( P 2 ) = - i , 

jr — m^ — ll(j? ) 
which defines the polarization function to be 

U(p 2 ) = (p 2 -ml)(l-f; 2 (p 2 )) . 

It has the necessary properties IT(p 2 = m 2 ) = 0, dU/dp 2 \ p 2 =m 2 = 1 — f~ 2 (ml) = because 
fdm 2 ,) = 1. Thus D R (p) can be viewed as the /u// propagator in this model. 

The WTI requires a relation between the full propagator and the irreducible vertex: 

(p " pO^,^) = e&r (^(p) - ^(PO) • 
Multiplying with the full propagator on each side gives the requirement 

(P - p')»D R (p')T» rr (p\p)D R (p) = eQ« (D R (p') - D R (p)) . (26) 

Requirement (^) can be satisfied with the following choice for the irreducible vertex: 

TUp',P) = (ml-A 2 ) 2 x 

r^(p',p)D c (p)DM + D n (p')T^(p',p)D n (p) + D„(p')D c (p')r^(p>,p)] , 

where T^(p',p) and T% (p',p) are the individual vertices of the free particles. This is 
pictorially represented in fig. |3] and can be shown as follows: 

(p-p') / ,D R (p')Vf rr (p',p)D R (p) 

= (ml - A 2 ) 2 eQ. [D c (p')(D-\p) - D- 1 (p'))D c (p)D^(p)D c (p) 

+ D c (p')D^p'){D-\p) - D-\p'))D^p)D c (p) 
+ D c (p')D n (p')D c (p')(D;\p) - D-\p'))D c (p) 

= (ml - A 2 ) 2 eQ n [D c (p')D„(p')D c (p') - D c (p)D ff (p)D c (p)] 

= e&(D R (p')-D R (p)) , 
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using the definition of D c (p) in (|24|) and the fact that the free vertex fulfills the WTI with 
the free propagator (pop. 

The naive picture, developed in the previous subsection, to couple a photon to each of the 
charged particles then very naturally emerges from the requirement of gauge invariance. This 
simplifies the numerical treatment of such processes considerably because the self energy of 
the pion need not be constructed explicitly; only free propagators occur, and, therefore, the 
pole structure is much more transparent. 

On the other hand, the "reduced" formalism provides a simple way to prove gauge 
invariance of the model. This is demonstrated for the exemplary case of a nucleon-pion- 
loop, where, in order to reduce the effort, the pseudoscalar coupling is used. The results 
also hold for the pseudovector case if the contact terms are taken into account. For the 
proof only two diagrams need to be considered: the photon couples to the internal nucleon 
(fig. 0b), and the photon couples to the internal pion (fig. |c). The pion propagator must 
be replaced by the reduced propagator and the vertex correspondingly by the irreducible 
vertex. For the moment the isospin factors are neglected, at the end they will be considered 
for protons and neutrons separately. With q fJL = (p — p')^ on gets for diagram [|c 

r d 4 k 
?/j J -7^T9l5iS (k)g-f 5 iD R (p' - k)(-i)T? rr (p' -k,p- k)iD R (p - k) 

using (26) 

= e J ^91,iS (k)g l5 iQ n (D R (p' - k) - D R (p - k)) . (27) 

(fZTj) is equal to the difference of self energy diagrams. But since Q n is still left in the 
expression, it is only that part of the self energy which is due to charged pions. The 
remaining contribution of the neutral pion comes from diagram [|b: 

r d 4 k 
% J 77^u9l5iS (p' - k)(-ieQ N Y)iS (p - k)g^\D R (k) 

= ~e j '-0^91 5 S o (p' - k)Q N (So\p -k)- S \p' - k))S (p - k)g l5 D R (k) 

f d 4 k 
= ej j^ri9j5iQN(So(p' -k)- S (p - k))g^iD R (k) . (28) 

15 



For an incoming proton only the 7r + n-loop contributes to fl27|) and gives an isospin factor 
2, the 7r°p-loop contibutes to (^) with an isospin factor 1. Adding both, one finds exactly 
e(£(j/) — S(p)). For an incoming neutron the 7r~p-loop contributes to both diagrams with 
a factor of 2, however, since the charge operators are present, there is a relative minus sign 
between ( |2"TD and (|2lf), such that they cancel exactly. The 7r°n-loop contributes to neither 
diagram. 

The WTI requires for the nucleon 

q,T% N ^p',p) = eQ N (S-\p) - S-\p')) 

= eQ N (^y ■ p-m- S(p) - 7 • p' + m + X(p')) 
= eQ N (j-q + Z(p')-Z(p)) . 

Obviously the direct term (fig. |2|a) accounts for the 7 ■ q whereas the vertex corrections 
exactly make up for the self energy. This proves gauge invariance up to the desired order 
in the strong coupling constant. The proof for pseudovector coupling or for the case of an 
internal A follows the same scheme as outlined above. 

With the presented method the electromagnetic vertex correction can be constructed 
from the self energy given to the desired order in the strong coupling constant. It is given 
by just adding external photon lines to each propagator corresponding to a charged particle. 
If a cutoff function is needed it can be absorbed in the "reduced" formalism, which is 
powerful enough to allow a proof of gauge invariance, on one hand, and simple to implement 
numerically, on the other hand. 

IV. ROLE OF VECTOR MESON DOMINANCE 

The idea of vector meson dominance was first introduced by Sakurai ||. It supposes 
that the photon couples to the hadron by first converting to a vector meson (a p-meson 
specifically), which then couples to the hadron. The idea was first investigated for the 
electromagnetic form factor of the pion, where it turned out to be successful. For time like 
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momentum transfers the structure of the p-meson propagator shows up clearly. Furthermore, 
based on the assumption that the p-meson is the "gauge-boson" of local isospin rotation, 
Sakurai introduced the vector meson universality hypothesis, which states that the p-meson 
couples to all hadrons with the same universal coupling constant g p , which is equal to the 
p-7 coupling constant g pl . If this was true, the electromagnetic form factor of the proton 
should have the same shape as the one of the pion. 

From electron scattering on nucleons one knows that the space like form factors are 
well described by the dipole fit [p]]. This is not in agreement with the simple assumption 
of VMD, which always results in monopole form factors. For this reason a number of 
resonances besides the p-meson are introduced. Their coupling constants to the nucleon 
are determined by a fit to the data |2|||, in a so-called pole fit. In the calculations of ref. 
0, furthermore, information of 7riV-scattering is included. By constraining the coupling 
constants it is possible to eliminate the leading monopole term and thus to obtain a dipole 
shape 

These pole fits can be analytically continued into the unphysical region (nucleon on 
mass shell, time like momentum transfer of the photon below the NN threshold). A rich, 
dominant pole structure shows up, which - if it were observable - would contribute a direct 
proof of VMD for the nucleon. The pole fits also predict the behaviour of the form factors in 
the time like region above iViV-threshold. Recent data taken at LEAR |4j] show disagreement 
with the fit obtained by || on the basis of the VMD hypothesis. 



It was pointed out by Holder et al. [15], that besides the poles in 7r + 7r~-annihilation, 
also the process NN — ► 7r + 7r~ is important to describe the form factor of the nucleon. The 
failiure of the pole fit of |J most probably is linked to this. It is thus an interesting question 
to investigate if 7riV-scattering influences the picture of VMD considerably. 

Since the information is hidden in the "unphysical region", the experiments must be 
performed using off-shell nucleons. It must thus be checked, if the off-shell dependence of 
the form factor interferes with the VMD picture. In this section a scheme is set up to model 
VMD for pions as well as for nucleons, which in section [V| will be combined with the off-shell 
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information. 

A. Coupling of Photons to Pions 

The electromagnetic form factor of the pion can be parametrized as 

2 
Tfl 

FM 2 ) = 2 ■ /, ■» j , (29) 

rrip — vm p i p{q z ) — q z 

where m p = 0.77 GeV is the mass of the p-meson, and 

T p (q 2 ) = 0.2458 (g2 ~ 4 2 m ' )3/2 e(g 2 > Ami) 

Q 

describes the decay of the p-meson into two pions, which is with more than 99 % the 

dominant decay-channel. 

The coupling that leads to such a form factor can be modeled by the following Lagrangian 

density. Since the only processes of interest for this discussion are those competing with the 

coupling of a photon to the hadrons, only the the neutral vector mesons p° and u will be 

considered. 

W + i(<W>" + e^)T 3 + g^) tt] 
_I m 2 - t ~ _ 1^/- + c pi + C 7 (30) 

4*, - ~G^CC + im^'* - ^-F^GZ , (32) 



4 M i/ u 2 w M 2 ^ u 



where G£„ = <9 M/ cv - cip M , G^ = c^a> v - d u u p , F pv = d^A u - d u A p , g pl and g pn are the 
coupling constants of the p-meson to the photon and the pion respectively, g^ and g uyn for 
the a»-meson accordingly. 



Note that with the tensor coupling in the Lagrangian (31, |32[) the photon-vector meson 
vertex turns out to be proportional to q 2 . In this case the contributions of the vector 
mesons to the photon polarization function vanish for q 2 = 0; the photon thus remains 



massless very naturally, whereas in the other case more involved constructions are needed 
to keep the photon mass zero |J . Also the coupling of the field tensors is gauge invariant 
by definition. 

For the form factor one gets 

2 2 

F w (q ) = 1 H : , 2 , 1: H j : , „. „- . (33) 

9 P ~i m P ~ irn p T p (q 2 ) - q l g wl m% - im u T w (q z ) - q 2 

Complete vector meson dominance (VMD) assumes that g p7T = g pl = g p || and neglects 

the coupling to the tu-meson completely. Under this assumption ([33]) reduces to ( ^9f) . The 

constants g p7T and g P7 , however, can be inferred from the decay-properties of the p-meson 



p —> 7r + 7r _ and p — >• 7 — ► e + e~. By converting the measured widths [29| into coupling 

constants one finds g p7T = 5.9 and g pl = 5.1. So the universality criterion is almost fulfilled, 

but not exactly. 

To calculate the ratio g un / ' g lx , 1 one uses 

9 (ml 9^ Tp^TTTr • T p ^ e + e - _ 1 
9l«/9l-y IV^™ • r^_ e +e- 0.017 2 

This shows that the u;-contribution to the form factor is much smaller than that of the 

p. Nevertheless, it shows up in the form factor because of the rather small cu-width which 

results in the structure on top of the wide bump stemming from the p-meson (see fig. f|). 

The picture that emerges is thus a little different from complete VMD: best agreement with 

experiment is reached if one relaxes the assumption of vector meson universality and allows 

to couple all particles, the photon and the vector mesons, to the hadrons with the appropriate 

coupling constants. In principle also the phases could be chosen to correctly describe the 

p-uj interference |1J . However, for the arguments in section the current agreement of data 

and theory is sufficient. 

B. Coupling of Photons to Baryons 

In this subsection the previously developed method will be carried over to baryons, and 
to the proton in particular. As stated in section |TJ the electromagnetic form factor has an 
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isoscalar and an isovector piece. By the VMD hypothesis these are related to the isoscalar 
cu-meson and to the isovector p-meson respectively. If one assumes a Lagrangian density 

C = $(j.(id- g pN T 3 p° - g uN u> - -(1 + r 3 )eA) -rn\v 
+£ + L F^F 



with £ p7 and £ W7 from (|3T|, |32|), one finds 



Ff P U tm (W, W; q 2 ) = Ffm\ W- q 2 ) ± F$(W, W; q 2 ) 



+F*':(W',W;q 



2N Q 



2 



ml - im^T^q 2 ) - q 2 



±F(*(W\W-q 2 ) 2 , f . (34) 

' p m 2 p - im p T p (q 2 ) - q 2 

F?J and F?y contain the direct coupling to the photon, F?J, contain the coupling to the 

vector mesons with the appropriate coupling constants. The structure of this Lagrangian 

density is guided by a gauge principle and thus minimal coupling for the vector meson fields. 

Additionally there exists a tensor coupling of the form 

- a^ v G p - a^ u G w 

C tensor = g pN ^K p —-^^ + g ujN ^^—-^^ . (35) 

4m 4m 

As a useful example the most simple case with the assumption of no further substructure 
of the nucleon will now be discussed. It is defined by the Lagrangian given above and neglects 
all couplings to further mesons like pions. Section [V| will combine the picture developed in 
section [Ifl| with the ideas given below to the complete scenario. 

For this simple example one finds Ff'g(W',W;q 2 ) = Ffy{W' , W; q 2 ) = 1/2, 
F?J(W', W; q 2 ) = g^N I g^ and F?J(W', W; q 2 ) = g P Nj ' g^- In principle there will be contri- 
butions to F2 by the tensor-coupling of the vector mesons, but these contributions cannot 
account for the anomalous magnetic moments of proton or neutron since they are weighted 
with q 2 and thus do not appear for real photons. So the anomalous magnetic moments are 
genuinely due to the inner structure of the nucleons. 

The ratios of the coupling constants can be determined with the help of meson-exchange- 
potentials for the nucleon-nucleon-interaction, if one takes e.g. |32| one finds g P N/ ' g^N ~ 0.26. 
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From the ratios of the decay widths one infers g n / 'g w1 = JT w ^ e + e - /T p _, e + e - ps 0.3. These 
ratios are remarkably close to each other, and will in this crude approximation assumed to 



be equal. Also from [32] one finds g P NJ 9(m ~ 0.47 which is surprisingly close to 1/2. It must 
be emphasized, even though this is only a qualitative discussion, the obtained numbers are 
close to what is obtained by other authors (table |). So as a summary one has for the proton 



2 1 (f 1 q z 

-f l,proton\Q ) J- + 77 n ■ r> T~2\ 2 



2 ml - im^T^q 2 ) - q 2 2 m 2 - im p T p (q 2 ) - q 2 

For space like q 2 and under the further assumption m p « m u « my one finds Fi proton (q 2 ) = 
my/(my — q 2 ), and F^ neutron (q 2 ) = 0. This establishes the commonly believed VMD hypoth- 
esis. The result falls short to explain the well established dipole fit to the electromagnetic 
form factor |jj since it is only of monopole structure. 

Concluding, one can state that the pion is essentially structureless, and that probably all 
of its size, if tested electromagnetically, is due to the p-meson. The nucleons obviously have 
structure besides the one due to VMD. This is very much in the spirit of Iachello, Jackson 



and Lande fL6 |, who introduce a further function to take into account the shortrange part 



of the interaction, or of the two-phase model of Brown, Rho and Weise JT7j , who explicitly 
introduce contributions of the quark-core. Also a lot of work has been done to model the 
structure of the nucleon by taking into account the meson cloud in terms of the exchange of 



virtual pions flT3|,|i"4" 



C. Gauge invariance and VMD 

At first sight VMD spoils gauge invariance, since even for the free vertex and propagator 
the WTI is no longer fulfilled in the simple model developed above: 

ml 

<^r M = e7 ■ q— + e7 -q . 

my — q z 



This can be restored by using the same technique as for the pion in section |HB|. Instead of 
the 7 M -coupling for the vector mesons a modification can be used: 
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Knn = 9vn ir - Q^j ■ (36) 

This vertex changes only the longitudinal part of the electromagnetic vertex (thus contributes 
only to F 3 ) since i ll T VNN = gvNJ^I 11 and restores the WTI since now q^XvNN = 0- 

V. RESULTS 



Using dispersion techniques, Hohler and Pietarinen [15] pointed out the relation between 
pion-scattering phase shifts and electromagnetic properties of nucleons. Because of a lack 
of higher lying resonances, it cannot be expected that a simple loop expansion like the one 
developed here can reproduce the required phase shifts for 7riV-scattering in the entire energy 
regime under consideration. However, from the Lagrangians ([H|) and (^) the phase shifts 
in the Pll and P33 channel can be successfully calculated close to threshold ||34j| . On the 
other hand, the mi — > 7T7r-amplitude in the J = T = 1-channel can be best described by the 
p-resonance. All important low energy thresholds are included for W > m by the inelastic 
channels in the loop expansion, as well as the thresholds for q 2 > 4m 2 by using the correct 
momentum dependent p-decay width as shown in fig. |5|. So it is natural, and no double 
counting is involved, to combine the loop expansion approach (section |ITl|) with the idea of 
vector meson dominance (section [PV|) . 

The form factors are thus constructed in three steps. First, all diagrams of fig. H are 



calculated with the Lagrangians from section [II. These diagrams form the contributions to 



Ffg(W', W; q 2 ) and F?y(W r , W; q 2 ). Next, these diagrams are calculated for a coupling to 
the neutral vector mesons instead, using the respective experimentally determined coupling 
constants for the p°- and cu-mesons to the hadrons and the photons, resulting in the con- 
tributions to Ffp(W, W; q 2 ) and Ffj(W', W; q 2 ). Combined as in fl34|), these contributions 
give the complete irreducible vertex. Solving a system of linear equations (see appendix) 
finally yields the irreducible form factors. The full form factors are obtained by applying 
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The relevant coupling constants are determined by experiment and symmetry consider- 
ations and are given in table [TI]. The ratios between the couplings of the pion to nucleons 
and deltas based on the SU(2)xSU(2)-considerations of |35| are carried over to the p- and 



u;-meson. The absolute values of the pNN and luNN vertices are taken from ||32|| . The 



AA/7 coupling constant G\ is chosen to lie in between the two values from | 2"q| . So, besides 
A, no free paramter is involved. Because the cutoff only regularizes the divergent term the 
results still need to be renormalized. In this paper the following renormalization scheme is 
employed (quantities with superscript R denote renormalized quantities): 

W* R (W) = W*{W) + c w , m* R (W)=m*(W) + c m , ff° R = Zff s . 

The numbers cw, c m and Z are constants, independent of q 2 or W; they can be determined 
from eq. fllTD- One finds 



/;:V;^»h r *- m " 



W — m 
For the renormalized theory one needs 



, p stands for proton . (37) 

W=m 



r*R ~~*R TT/*^ I iv^*^ 



s-\p)s ( P ) = A + ( P r Tir m +A-( P r 



W — m ~ W + m 

u~l 1 = A+ + A ~ • (38) 

Comparing the coefficients of the projection operators, it turns out, that in this limit 
w w Zm — goes to 1, which is equivalent to the statement that the full propagator has a 
pole with unit residue at W = m. Therefore, one has W* (m) = m* R {m). It cannot be 
deduced, however, that W* (m) = m and m* R (m) = m, because for W — » m the projector 
A~ becomes zero itself, and thus there is no constraint on W* (m) + m* (m). 

For a pole of first order one has Res(f(z); z ) = lim 2 ^ 20 (z — z )f(z). With help of this, 
one can read off from (|37| ) that V/iV (m; 0) is the value of the residue of the propagator and 
must thus be equal to 1. This residue is adjusted with the wave function renormalization 
constant Z. Therefore, one must choose Z = l//^ p h (m; 0). It turns out, that the full 
form factors are renormalized automatically since, because of eq. (|19D , the wave function 
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renormalization constant drops out. In the actual calculation Z = 0.37 is found. For the 
self energies one gets W*{m) = m*(m) = 0.9m, thus £(ra) « 0.1GeV(7p — m), which is a 
reasonable value. 

A. On-shell Form Factors 

Fig. P shows the results for the full on-shell form factors F^iq 2 ) f° r proton (a,c) and 
neutron (b,d). Note, that in this case the full form factor is equal to the irreducible one. 
The squares are calculated from the dipole fit to Ge(q 2 ) and Gm{q 2 ), which describe the 
data for these momentum transfers up to a few percent in this energy range Jl|. In (b), in 
addition to the dipole fit, information about the deuteron form factor has been taken into 



account. The two curves of squares correspond to the extreme parametrizations of |36| to 
indicate the theoretical ambiguities in the description of the NN-potential for the deuteron. 
A severe test of the model, F^ eutron (q 2 = 0) = 0, is well fulfilled. For the proton, 
furthermore, -F\ ++ is reproduced very well. The radius of the proton comes out to be 

dF+ + (q 2 ) 



(r 2 ) = 6 



=> (r 2 ) 1 / 2 = 0.81 fm . 

g 2 =0 



dq 2 
Fi~ + for the neutron is comparable with the data; it falls in between the theoretical uncer- 



tainties of the experimental analysis. As discussed in the toy model of section |TV| , the VMD 
contribution to the form factor accounts for a monopole shape only. Fig. ^ indicates, how- 
ever, that in the full model a dipole shape is obtained. Therefore, it is possible to conclude 
that in addition to the explicit treatment of the vector meson pole terms, the 7r-loop cor- 
rections are thus necessary to obtain the experimentally observed dipole shape of the form 



in a 



factors. This conclusion has been discussed previously by Gari and Krumpelmann p3 
qualitative way, and is in agreement with the findings of early dispersion theory treatment 
of the form factors. These results give some confidence in the validity of the off-shell results. 
Fig. |5| c,d displays the behavior of F^~ + (q 2 ). The shape of the form factor is in agree- 
ment with the data in the momentum range plotted; however, it does not show the correct 
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behaviour at larger momenta. Furthermore, the model cannot account for the full size 
of the anomalous magnetic moment. The values found in this calculation are k v = 1.45, 
K n = —1.65, to be compared with the experimental values of Kp Xp = 1.79, K e * p = —1.91. As 



mentioned in section [IV B|, in the Lagrangian given in fl3~l"|) and fl32|), there is no contribu- 



tion of the p-tensor coupling to F 2 for real photons. This is in contrast to |14}|, where, by 



multiplying the bare vertex with k p /2 x my/(rriy — q 2 ), an anomalous magnetic moment is 
induced. Such a procedure is highly questionable because it can be argued that the tensor 



coupling of the p- meson to the nucleon is just due to the loop corrections discussed here |]37 
To avoid any double counting the tensor coupling must, therefore, not occur in the direct 
diagram. In the present calculation, it has been taken into account only for the diagrams 
with internal radiation to simulate the higher order corrections in a schematic way. This is a 



crude approximation of the iteration scheme devised in J38] . The tensor coupling influences 
only the g 2 -dependence of F 2 , not its magnitude at q 2 = 0. A number comparable to the 
magnetic moment calculated here can be obtained by dividing the result of fLM with meson 



cloud by their result for the bare vertex, giving roughly 2/3. Also ITS] gets a number which 
is too small (k p ~ 0.5). 

Fig. [F] shows the irreducible form factors for projections to negative energy states. f± 
for the proton is not constrained to 1 by (|38|) . Its value of 0.9 at q 2 = is directly related 
to the scalar and vector self energies. Neither is fi for the neutron constrained to 0. The 
magnetic form factors show a change of sign as one goes to higher q 2 . Their values at q 2 = 
are drastically different from f2 + {q 2 = 0). The full form factors F^ can be obtained by 
rescaling /^ by 1/0.9. Even though these form factors never play a role in experiments, 
they already indicate that all theoretical calculations that rely on the asumption F^ + =F^ 
are incorrect. 



The models |]13| , |T4l , |2"5| all approximately agree with each other for the electric form factor. 
Since their magnetic moments are off by more than 50 % from experiment, an important 
process must have been missed in all these calculations. There exists strong experimental 
evidence, that a large part of the magnetic moment is due to spin flip transitions, especially 
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to the A- resonance fL8[ . The present calculation shows, that about half of the magnetic 
moment is carried by the diagram with the A — > Nj decay mode, which is dominantly Ml, 



a mode that was not included in |TJ^[rj|,[2IJ . 

It must be stressed, that the underestimation of k does not depend on the cutoff A used 
in the calculation. Fig. ^ displays the A-dependence of k p and K n . The magnetic moments 
are optimal in the 1 GeV region, giving confidence in the combination of the employed 
regularization and renormalization scheme. One clearly finds for the larger values of A 
that the magnetic moment decreases. Also for small A there is the onset of a decrease, 
which shows that the choice of the cutoff parameters cannot solve the disagreement. This 



is in strong contrast to |37|, where the tensor coupling increases monotonically with A. 



The difference can be traced to the renormalization procedure of [p7|| , where a subtractive 
renormalization is used, which is valid only, if the wave function renormalization constant 
is close to 1; this certainly is not the case in the present calculation [Z = 0.37). Note that 



the agreement of state-of-the-art soliton models that include the A-resonance [12|] with the 
magnetic moment is of the same quality as the agreement of the model presented here. 
Bos et al. P5| dress the nucleon with a scalar /isoscalar cloud. They use a cr-meson with 



a mass of m a = 0.8 GeV. Since in this case the photon only couples to the nucleon, this 
gives a large weight to the core contributions. Still, the magnetic moment is only around 
0.7. If a contribution like this is included in the present calculation, there is some effect, 
which is displayed by the dashed line in fig. §. The parameters chosen here are m a = 0.56 
GeV, g a = 11. The effect is due to a reduction of the wave function renormalization Z. This 
increases the magnetic moments somewhat, however, without coming anywhere close to the 
experimental point. 

B. Space Like half off-shell Form Factors 



Figs. and ITU] show the off-shell dependence of F x + (W; q) and F£ (W; q ) for the 



proton. In order to see the effects of the pure loop corrections all contributions due to VMD 
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were switched off in the upper panels. Part of the purpose of the following discussion is to 
show, that not only in the time like electromagnetic form factor one can see the influence of 



the NN — > 7T7r-scattering |15[], but also in the half off-shell form factors. 

The left column of fig. |9] shows -Fi ++ . For W < m + 2m,,- the slope of the real part of F^ + 
increases with W, corresponding to an increasing charge radius, which is in agreement with 
|13| , [T^| . Above m + 2m n the radius stays more or less constant, the real part falls rather 
linearly with — q 2 . The imaginary part reflects the possibility that the incoming nucleon 
is far enough off-shell to decay into final states of one or two pions plus a nucleon or a A 
during the electromagnetic scattering process. Because of the various inelastic thresholds, 
the imaginary part is expected to depend on the phase space of the decay products. This 
phase space dependence, however, is mixed with self energy corrections due to the WTI, 
which themselves have imaginary parts. Thus the imaginary part of F^ + does not show a 
clear behaviour as q 2 and W increase. The absolute value of the form factor including VMD 
shows only little dependence on the off-shellness of the incoming proton for F±~ + . 

For F^~ (right column of fig. |9|) the situation is different; above the 27r-threshold it 
rises fast. Even after including VMD, the changes due to off-shellness are clearly visible. 
Calculations of cross sections for electromagnetic processes including off-shell nucleons and 
the correct form factors must be performed in order to show the importance of the negative 
energy components. 

For the magnetic form factor -F 2 +± (fig. |T0|) there is a dramatic change in shape above 
the 27r-threshold in the real part, whereas the imaginary part mainly seems to grow with 
phase space. At q 2 = the real part of -F 2 +± develops a maximum as a function of W at the 
27r threshold; the absolute value, on the other hand, increases rather smoothly. 

A detailed analysis of the contributions of specific diagrams to the full form factor reveals 
that the above effects are really caused by the N* — > ttttN channel, where N* denotes an 
excited nucleon. Fig. |11] shows the decomposition for W = 0.939 GeV and W = 1.4 GeV in 
comparison. The full symbols show contributions of radiating pions (fig. 0c-e), i. e. the meson 
cloud, the open symbols correspond to radiating baryons (fig. §b), i. e. the core. Squares 
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describe diagrams with an internal nucleon, circles represent diagrams with a propagating 
A. The solid line is the contribution of the direct diagram (fig. |2]a) and is counted as a core 
contribution. Contributions from the diagram of fig. |2|f are indicated by dashed lines. 

It is quite important to notice that for iq for both values of W the cloud contributions 
appear with different sign, leaving a more dominant core, whereas for F 2 the core contri- 
butions, at least for W = m, tend to cancel. So while iq is core dominated, F 2 is more 
influenced by the meson cloud. The most important contribution to F 2 is the transition of 



the A, in agreement with the analysis of the Gerasimov-Drell-Hearn sum rule fig]. 

An obvious change in iq ++ , if one goes off-shell, is the increased importance of the meson 
cloud from 15% to 35%. This is clearly an important sign that the mesonic excitations must 
be carefully treated in all models for off-shell form factors. Already in the on-shell case the 
contributions with a propagating A are of some importance. 

From the slopes one can deduce that the contributions of radiating pions (cloud) reach 
further out than the ones of radiating baryons (core). Since the different contributions show 
different slopes at q 2 = 0, which is especially true for the off-shell case, it is not possible 
to give an adequate parametrization of the q 2 dependence of the form factor using a single 
parameter, if one wants to maintain a cloud/core picture to describe the extension of the 
nucleon, like in the semiphenomenological models of |16|JT7|] . This is even more clearly visible 
for the magnetic form factor. 

In going from the on-shell point to 1.4 GeV, the major changes come from the diagrams 
^|c-e. For iq ++ these fall off faster for 1.4 GeV; for F 2 ++ the shape is completely determined 
by these diagrams, all other contributions stay rather constant. Since these diagrams are 
the only ones that allow for the 27r-decay, it must be concluded that the N* — > tcttN channel 
is responsible for these effects rather than the N* —>■ 7rA channel which is already open, too. 
This also indicates that the effects described above are not spurious and will remain if the 



A is treated in a more realistic way by allowing for a finite decay width (cf. section [Til A|) 



2N 



C. Time Like Form Factors 



Figs. 12 and 13 show the electric and magnetic form factors for the proton, respectively, 
in the entire range covered by this calculation. The contour lines show \og\F^(W 2 ;q 2 )\. 
Various thresholds are indicated by dotted lines. The region, which is not accessible by 
experiments ('unphysical region') is indicated by the gray area; all dispersion relation ap- 
proaches rely on analytical continuation into this region. 

It must be emphasized, that at this level all inelasticities up to the 2tt channel are 
included in the vertex either by the p — > tttt decay, or by N* — ■> X + mc (n = 1, 2). While 
F\ shows almost no changes as one goes off shell, for F 2 the influence of inelastic channels 
remains visible. Especially the onset of 2n production influences the form factor. 

It is due to the WTI that F\ does not change very much when going off-shell. Each 
off-shell effect in the vertex is nearly canceled by the self energy corrections (|I9|) , which have 
the same thresholds. Interestingly, this cancellation seems to be independent of q 2 . This is 
not entirely true (see last subsection), however, the q 2 dependence is so much dominated by 
VMD that other effects, which depend on q 2 , are hardly visible. This is a very important 
result. If VMD exists for nucleons as well as for pions, it will clearly be observable, since it 
is not reduced for off-shell nucleons in the case of Fi, and is still dominant for F 2 . 

For W = 2 GeV fig. [13] displays the scalar and vector contributions of the direct coupling 
graphs and of vector meson graphs to the form factors. From this information the proton 
or neutron form factor can be constructed. For F\s/v it shows a smooth behaviour; the 
minimum in the vector imaginary part is due to the 7r + 7r~-channel. The pure photon channels 
are clearly least important as compared to the vector meson channels. The figure shows the 
broad rho- and the narrow omega-resonance. Since both mesons couple to the nucleon 
with about equal strength, the uj contribution is an order of magnitude larger than the p 
contribution at peak level because of the cj's small width. Note that in the case of F 2 for 
the vector meson contributions the role of imaginary and real part has changed, indicating 
that the imaginary part of the loop expansion without VMD is significantly larger than its 
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real part. This is consistent with the situation in the space like sector (fig. pi]) and gives rise 
to interference. For F 2 the contribution of the p-meson if of larger importance than for F\. 
To show once again that VMD can be observed even in an experiment involving a half 



off-shell vertex, a cut of fig. [12] is plotted in fig. |15|. The absolute magnitude of the single 
contributions in the photon-, to- and p-channels as well as the absolute magnitude of the 
coherent sum is displayed for time like momentum transfers for the proton and the neutron 
at W = 2 GeV. Despite the small width of the tu-resonance it is not possible to resolve the 
p-contribution in F\. The sum is almost exclusively exhausted by the ^-contribution, while 
the p on the other hand is only visible as a broad background. On the other hand, for F 2 
the p-component is not much smaller than the tu-component, therefore, subtle interference 
effects show up between the p- and u;-channel, which might be accessible in experiments on 
both protons and neutrons. 

Experiments on dilepton production in p + A reactions, where the 7r + 7r~-annihilation 
channel is surpressed, may thus offer a chance to study medium modifications of the cu-meson 
from the electromagnetic properties of the nucleon in the same way as medium modifications 
of the p-meson from electromagnetic properties of the pion. Because the nucleon's form 
factors are larger than the one of the pion (at peak position \F K \ « 7), it is visible even in 
heavy ion collisions above the background from -k + tt~ annihilation ||; however, a very good 
resolution of the experimental apparatus is required to resolve it. 

VI. SUMMARY 

To study the production of dileptons in high energy nuclear reactions one needs infor- 
mation about the half off-shell time like electromagnetic form factors, which are until now 
unknown. We have, therefore, constructed a dynamical model based on a hadronic frame- 
work to calculate the electromagnetic form factors for momentum transfers of the photon of 
-1 GeV 2 < q 2 < 1 GeV 2 and for nucleons with -1 GeV 2 < p 2 < 4 GeV 2 . 

Starting from an expansion of the nucleon propagator in pion loops, the electromagnetic 
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vertex is constructed by inserting external photons to each charged particle line, thus obeying 
the constraints due to the WTI. Together with the concept of VMD, this approach includes 
at least schematically all findings of the spectral analysis of the late 1960 flT5| . It is possible 
to maintain a cloud/core picture that proved successful in semiphenomenological models 
T5| , [TT[ ] . It is important to include the A-resonance to reproduce the magnetic properties, 



which is in accord with an analysis of the Gerasimov-Drell-Hearn sum rule [18]. The coupling 
constants and the cutoff are chosen in agreement with meson exchange models and symmetry 
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The momentum dependence of the space like on-shell form factors is reproduced. The 
charge radius of the proton is found to be r = 0.81 fm. The electric form factor of the neutron 
falls in between the uncertainties of the data analysis. The magnetic moments are k p = 
1.45 and K n = —1.65, which is much closer to the experimental value than in comparable 



calculations [p^ , ^4j] , but still too small. The better agreement can be traced to the A — > N'j 
decay process, which occurs in neither of the above cited works. A decomposition of the 
form factors shows that for F\ the core contributions dominate, while for F 2 the cloud is 
more important. It also shows that inelastic thresholds influence the form factors. However, 
for F\ the dependence is only weak because it is compensated by self energy corrections 
required by gauge invariance. For F 2 threshold effects remains visible. 

The weak off-shell effects enable one to study VMD in the experimentally accessible 
region of the (W, q 2 ) plane. It turns out, that only the signal of the w-meson can clearly 
be extracted, the p-meson contributions only result in a broad background. As an effect 
one has the possibility to study medium-effects on the u very clearly by measuring the 
electromagnetic form factor in p + A- and even in heavy ion collisions. 

It is found, that calculations which are based on the assumption that F^ + = F^~ must 
be rechecked. 

The model for the vertex used here certainly leaves room for improvements. Taking into 
account higher resonances, especially of spin 3/2, as well as heavier mesons will help to im- 
prove the magnetic moments. The finite decay width of the A also needs to be incorporated. 
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But since the most prominent threshold effect is due to N* — > Nmr, it is expected, that the 
off-shell behaviour of the form factors will not change appreciably. 
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APPENDIX A: NUMERICAL DETAILS 

The sum of all Feynman diagrams in fig. ^] results in the irreducible vertex. To obtain 
the form factors from the most general vertex (|J) the projection on positive /negative energy 
must be performed for in- and outgoing nucleons: 

A'V)rV,P)A 8 (p) = A s '(p>)J2f? s (W',W;q 2 )0!tA s (p) . (Al) 

t=i 

Performing traces over contractions of the vertex with three linear independent 4-vectors 

Vy = P M = p /M + p M , V2=(f = p" 1 — p M , U3 = 7 M gives 

f t = Tr v?A s '(pi)r,A s (p) = J2 a l3 F°' s . (A2) 

i=i 

The same procedure is performed with the sum S^ of all contributing Feynman diagrams: 

T, = Tr < A*V)^A s (p) . (A3) 

To obtain the form factors one must equate ( |A2|) and 



T i = f i = Y g a ii F^ ■ (A4) 

i=i 

While the expressions Tj can be obtained analytically, the quantities T, must be calculated 

numerically. They can be decomposed into integrals of the following form: 

d 4 k 



Ti = T 



E 4%w\p 2 , k 2 )(p f knpknk 2 r/N d (p\p, k) 



The index d labels the different diagrams contributing to 5* M . Nd(p',p,k) contains the de- 
nominator of the propagators in diagram d. 

The coefficients c^f' v s w {p' ,p 2 ,q 2 ) as well as the o^ in eq. ( |A"2]) are calculated using the 
high energy package of REDUCE |39| . 

All integrals are of the form 

d A k {p'k) u {pk) v {k 2 ) w 



I{p',p) 



(2tt) 4 jv 
v 



II ((a p - k) 2 -m 2 + ie) 
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where a p and m p are the specific momenta and masses of the propagators in the diagrams. 
The number N of factors in the denominator depends on the choice of monopole or dipole 
cutoff and ranges from 4 to 7. The integrals are performed numerically. 

I(p',p) is solved in the restframe of the outgoing nucleon, which always exists in case of 
half off-shell kinematics. The integration of k° is done by contour integration around the 
poles of Il^Li(( a p — k) 2 — m 2 + ie) _1 . Let kf denote the poles in the upper half plane. Since 
not all poles are of first order, for each k the residue must be calculated by taking numerical 
derivatives: 

r d 3 k n 

^.^feEH ) , (A5) 



1 cf 1 * -1 



{p'k) u {pk) v {k 2 ) w {k°-k°) 

N " 

I! ((Op- k) 2 -ml + ie 
p=i v 



fc°=fc? 

where n-i is the order of the pole at fc°. For space like momentum transfer the sum of the 
residues is a well behaved function of k. However, in the time like region it still has poles 
because of the physical inelasticities for q 2 > 4m 2 and p 2 > (m + m^) 2 . These poles are 
treated with a subtraction technique, which will be described below. Since the momenta 
p' and p can be chosen such that / is invariant under rotations about the z-axis, the 
integration is trivial. The remaining integration is two dimensional. From now on k denotes 
\k\. 

//* f(k t] 
dx / dk—— 1 — - with x = cos9 
_, { 9{k,x) 

Let ki(x) be defined by g(ki(x),x) = (i = 1,2). The ki(x) are complex functions of x. 
For x > xq the imaginary part becomes ±e. Below xq an ordinary integration can be used. 
Above Xq the integral I(p',p) splits into a principle value integral and an imaginary part. 



n g(k,x) 7 g R (k,x)Ui=i[h(x,e = 0) -k + ieRe(dki/de)] 
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p hlk f{k ' x 



g R (k,x) [kx(x,e = 0) - k][k 2 (x,e = 0) - k] 



_ i7i f(h(x),x) dk t /de 
g R {k 1 (x),x)k 2 {x)-ki(x) 
. f{k 2 {x),x) dk 2 /de 

+ » ^M»),«)M«)-W 9(tl( '" (A6) 

Since P J °° l/(k 2 — kfydk = 0, the principle value integral is treated as follows: 



oc 

p fdh 



g(k,x) 



OO 

dk 

o 



f(k,x) *f{ki{x),x) k 2 -k 2 {x) 



Q{ki(x)) 



g(k,x) ~{k 2 -kf(x) g{k,x) 

The integrand is now finite for all k. 

If p 2 > (m + 2m n ) 2 , then in some diagrams both functions ki(x) contribute to the 



imaginary part. In this case there is a remaining singularity of the type 1/yl — (x/xq) 2 at 
x = xq. The integral is thus converging, but numerically unstable. It can be treated by a 
trick similar to the one used above: 

^ dx= fvW-^dx + Mxo) 



2 J * I rt*4 — o™* 

\xq\ V ° ~\ x o\ 



Jjf) X |_ i \/ Xq X 



After having calculated all integrals in Tj, the equations ( |A4| ) are solved for the invariant 
form factors. 

Numerical inaccuracies can occur at various points. For example, in ( |A5|) only the sum 
of all residues falls off fast enough in k so that the integral converges. There is a delicate 
cancellation of the summands which can be shown analytically that is hard to reproduce by 
taking numerical derivatives. Furthermore, for large k the poles come near to each other, 
sometimes too close to take the numerical derivative with good enough accuracy. This is 
worked around by not evaluating the integrand at too large values of k. A cutoff is chosen 
dynamically if the integrand is small enough. For higher values of k asymptotic behaviour 
is assumed and the integral is solved analytically. Unfortunately under this assumption the 
WTI suffer. They are only fulfilled at the 1% level. 
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In (|A6|) problems arise when \xo\ is close to 1. Then too few gridpoints of the integration 
mesh contribute to the imaginary part. In these cases additional grid points are created and 
integration weights are redistributed such that the integral is performed over at least 10 grid 
points in cos ^-direction. 

The numerical difficulties are best under control for a very simple integration with equally 
spaced grid points of equal weights. To obtain an acceptable accuracy up to 400 points in 
radial direction and 100 points in cos 8 direction are used. 
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APPENDIX B: FIGURE CAPTIONS 

Fig.l: Expansion of the full nucleon propagator in first order in pion lines. The nucleon 
propagator is displayed by the solid single line, the A-propagator by the solid 
double line, the full propagator by the solid line with the fat dot, the pion 
propagator by the dashed line. 

Fig. 2: Loopwise expansion of the irreducible nucleon vertex, (a) is the free vertex, (b) 
and (c) come from coupling of photons to charged hadron lines of the propagator 
in fig. |l|, (d) and (e) arise from contact terms due to pseudo-vector coupling of 
pions to nucleon and A, (f) is the contribution of the decay A — > N'-f. 

Fig. 3: Additional diagrams to satisfy WTI if a cutoff is introduced in the nN- or the 
7r A- vertex. The double dashed line represents the 'propagator of the cop'. 

Fig. 4: Possible cuts of the loop diagrams. The cut in (a) corresponds to the decay of 
the off-shell internal nucleon to an on-shell nucleon and a pion, the cut in (b) 
corresponds to 7r + 7r~-annihilation if the photon has high enough invariant mass, 
or to two pion production if the incoming nucleon is far enough off-shell. 



Fig. 5: Electromagnetic form factor of the pion. The data are taken from [p3| , |3"0| , |3T 
the solid line is obtained using fl3"3] ) and adding a similar contribution for the $- 
meson. 

Fig. 6: F^(m,m;q 2 ) for (a,c) proton and (b,d) neutron. The solid line is the model 
calculation, the symbols represent experimental results as explained in the text. 
For F 2 the experiment is rescaled to the anomalous magnetic moment of this 
calculation. 

Fig. 7: f\2( m i m 'i<l 2 ) for ( a ; c ) proton and (b,d) neutron. 
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Fig. 8: The anomalous magnetic moment for proton and neutron as function of the 
cutoff parameter A. The solid and dashed lines display the result with and 
withot contributions of the a-meson respectively. 

Fig. 9: Half off-shell form factors of the proton in the space like region. The left panel 
displays F^ + (W; q 2 ) for several incident invariant masses. The curves with sym- 
bols are below the 27r-threshold, the curves without symbols are above. The 
two upper plots show real and imaginary part of the form factor without VMD 
contributions, while the lower plot shows the absolute value including VMD. The 
right panel displays F{ + {W;q 2 ) in the same way. 

Fig.10 : Same as fig. | but for F^ + (W; q 2 ). 

Fig. 11: Detailed analysis of the contribution of the single diagrams in the W = 1.4 
GeV case (right panel) compared to the on-shell case (left panel) for the proton. 
Lines with full symbols correspond to meson cloud contributions, lines with open 
symbols and the full line are core contributions. The dashed line represents the 
A —>■ N'-f decay contribution. 

Fig. 12: -F 1 ++ for the proton in the full range of applicability of the model. The contour 
lines are steps of 0.1 in logio\F\. Various thresholds are indicated by dashed 
lines. The experimentally unaccessible region is marked as the grey area. 



Fig. 13: Same as fig. [12], but for F s 



Fig. 14: Real and imaginary part of scalar and vector contributions to the electric 
form factor for W = 2 GeV. The full lines display the isoscalar contributions, 
the dashed lines display the isovector ones. Symbols indicate imaginary parts. 
The upper plots show the scalar and vector contributions of the direct photon 
coupling. The plots below show contributions of u- and p-mesons seperately. 
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Fig. 15: Decomposition of l-Fj^j 1 "! f° r the proton (left) and neutron (right) at W = 2 
GeV into components coming from the direct photon vertex, the p-meson vertex 
and the tu-meson vertex. The sum is the result of a coherent superposition of 
the single contributions. 
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TABLES 





9uin/9w^i 


gpN/gp-f 




Hohler et al. § 


0.52 


0.98 


includes ttN phase shifts in p already 


Dubnicka § 


0.376 


0.418 


remaining strength in higher poles 


Gari et al. [33 


0.377 


0.411 


remaining strength in direct coupling 



TABLE I. Comparison of coupling constants in VMD like polefits 



9NNit 
9NNlu 
dNNp 



13.45 , g NAn = 22.86 , ^aatt 
8.34 , g NAu> = 14.18 , g AAuJ 
7.26 , g NAp = 12.35 , g AAp 
G x = 2.5 GeV" 1 

Kp = , K u = U 

A = 1.2 GeV 



10.76 
6.67 
5.81 



TABLE II. Coupling constants 
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